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We consider the inflation produced by two coupled fluids in the flat Friedmann-Robertson-Walker 
universe. Different cosmological models for describing inflation by use of an inhomogeneous equa¬ 
tion of state for the fluid are investigated. The gravitational equations for energy and matter are 
solved, and analytic representations for the Hubble parameter and the energy density are obtained. 
Corrections in the energy density for matter inducing the inflation and the coupling with energy are 
discussed. We analyze the description of inflation induced by non-constant equation-of-state param¬ 
eters from fluid viscosity. The correspondence between the spectral index and the tensor-to-scalar 
ratio recently observed by the Planck satellite is considered. 
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I. INTRODUCTION 

According to recent observational data the current universe is accelerating [l| . There exists an early-time accelerated 
epoch after the Big Bang, called the ’’hot scenario” or inflationary period (for reviews see Refs. [3,3). During inflation 
the total energy and the scale factor both increase exponentially Q. This phenomenon suggests that the cosmic 
fluid has properties quite different from standard matter and radiation. The simplest description of the accelerated 
expansion of the universe is given by the dark fluid model (the Friedmann-Robertson-Walker (FRW) model), in which 
the cosmic fluid is taken to satisfy a homogeneous or alternatively an inhomogeneous equation of state. 

Various problems associated with viscous cosmology have been discussed extensively in Refs. dill. References 
[ 234 ^ were devoted to inhomogeneous fluids, of which a time-dependent equation of state (EoS) in the presence of 
viscosity represents a sub-class. The behavior of a non-classical inflation was considered in Ref. (26|. Inhomogeneous 
fluid cosmology may also be understood as modified gravity (for reviews, see Refs. 1^ 1^ 1. owing to the fact that 
it may be modeled as a gravitational fluid with an inhomogeneous equation of state [29j. Then, it is to be expected 
that following the unification of inflation with dark energy in modified gravity [30l| . one can achieve this unification 
in viscous cosmology also, including the presence of dark matter. 

Some examples of inhomogeneous viscous coupled fluids were considered in Refs. [sTI - ls^l . In the papers [H, 
a bounce cosmology induced by an inhomogeneous viscous fluid in FRW space-time was investigated. In Ref. 
different kinds of inhomogeneous viscous fluids were analyzed concerning the possibility to reproduce the current 
cosmic acceleration in flat FRW space-time and the presence of finite-future time singularities, while in Ref. [s^ the 
noticeable fact was proven that a viscous fluid with an inhomogeneous equation of state is able to produce a Little 
Rip cosmology as a purely viscosity effect. Various cosmological models of the coupling between energy and matter 
were considered in the review [s^. 

In the present article we will investigate cosmological models in the presence of inflation. Some variants of inhomo¬ 
geneous viscous coupled fluids, of interest in the inflationary regime, are considered. The influence of the interaction 
between energy and matter on the description of inflationary cosmology is analyzed. The slow roll parameters, the 
spectral index, and the tensor-to-scalar ratio are calculated. Restrictions on the thermodynamic parameters, needed 
to satisfy the Plank data, are obtained. The agreement between theoretical inflationary models and observational 
data is discussed. 
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II. VISCOUS FLUID MODELS FOR INFLATION 

In this section we will consider the early-time universe, applying the formalism of inhomogeneous viscous fluid 
in flat FRW space-time. We will describe the inflation in terms of the equation-of-state parameters, and the bulk 
viscosity. We will pay attention to the behavior of the Hubble parameter and the energy density in the beginning and 
in the end of the inflation. 

We will start from a universe filled with two coupled fluids: energy and matter in a spatially flat metric with scale 
factor a. The background equations are (see e.g.. Ref. [s^): 

( p + 3H{p + p) = -Q, 

/ Pi + 5H{pi + pi) = Q, (1) 

i H = —^{p + p + Pi + pi), 

where H = a/a is the Hubble parameter and = 8ttG, G being Newton’s gravitational constant. Further, p, p and 
Pi , Pi are the pressure and the energy density of dark energy and dark matter correspondingly, and Q is a function 
that accounts for the energy exchange between the fluids. A dot denotes derivative with respect to cosmic time t. 
The cosmological constant A is set equal to zero. 

As is known, there exists a number of studies on cosmology with multi-component viscous fluids, especially in order 
to explain the late-time cosmic acceleration - cf., for instance, Refs. [Tsl - I^ Is^. IdOl - ld^ . It is then most natural to 
allow for interaction terms between the components. This means physically that each component is to be regarded 
as a non-closed physical system, corresponding to an energy-momentum tensor that is not divergence-free in general. 
The requirement of vanishing four-divergence can be made only on the total energy-momentum tensor i.e., 

^Total- 1 / “ element in our analysis is to adopt a two-component model also in the inflationary region. One 

may ask: is there a physical reason for adopting this model in the very special inflationary epoch? To our knowledge, 
there is at present no basic reason. Rather, the model is introduced on phenomenological grounds, to make the theory 
also in the inflationary epoch similar to that commonly accepted in other epochs. It should here be emphasized that 
the state of acceleration is common for inflation and for the late universe, especially when the Big Rip is approached. 
In a future work we plan to extend the present analysis so as to unify inflation with dark energy, such that one 
component describes inflation where the other component effectively describes the dark energy epoch. In this way we 
hope to present a unified universe evolution with the use of a two-component fluid. 

Let us now write down Friedmann’s equation for the Hubble rate 

= -^{P + Pi)- (2) 

We will assume that there is no relativistic matter {w = 0), and we will assume, as usual, that the dark matter is 
dust. Thus, Pi = 0. We also put c = 1, as usual. 

The gravitational equation for matter reduces to 


Pi -|- iHPi — Q. 


(3) 


We will now study some cases of inhomogeneous viscous coupled fluids in the inflationary period. 


A. Fluid model with oj = —p/{p + p«), and viscosity proportional to H 


Let us choose the viscosity to have the form [44 1 


C(H) = /(iL), (4) 

where the star subscript refers to the end of the inflationary period. The function f{H) will in the following be 
adjusted to have a physically reasonable form. Correspondingly, refers to the end of inflation. The 

function ([5]) is assumed to vary slowly with time, as inflation can be realized when the viscosity is slight. 

Let us take the equation of state to have the following inhomogeneous form 


p = oj{p)p + C{H), (5) 

where the thermodynamic parameter is uj{p) = —pj{p -I- p*) [4^, and p* is the energy density at the end of inflation 

(p. ^ 
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First of all, we will find the solution of the gravitational equation for the energy component 

p + mp[l + w(p)] - 3H^C{H) = -Q. 

Let us suppose that the ratio r = pi/p is a constant, and rewrite Friedmann’s equation ([2]) in the form 


( 6 ) 


fc2(l + r)’ 

We will consider the following coupling between the fluids: 

Q = 9SH^, 


(7) 

( 8 ) 


where i5 is a positive dimensional constant. In geometric units, its dimension is [i5] = cm“^. 

Since there exists no fundamental theory specifying the functional form of the coupling between the fluids, the 
present coupling model has necessarily to be phenomenological [s^. 

Further, we will analyze the model in the asymptotic limit corresponding to the initial stage of inflation. 




(9) 


Then <c 1) = f{H), and the function f{H) can be taken to be proportional to H, 


f{H) = 0H, 


( 10 ) 


with 0 = 3/(fc^(l + r)). 

Taking into account Eqs. ([US, 10) the continuity equation ([6|) for energy simplifies to 


-H + H^ -^ 

3 0H^ 




= 0 . 




( 11 ) 


As we are considering the initial stage of inflation, we may put the density p* in the form p* = 0H?^ with Hm = H (tin), 
where tin is the initial time. In this approximation the solution of Eq. (El) becomes 


H(t) = (V^r^ + 1 - T)Fti„, 


( 12 ) 


where r = (3/4)iLin(t — tin). The Hubble parameter is positive, so we have an expanding universe. 

Now, we solve the gravitational equation of motion ([3]) for matter with the coupling terms ([5]) and the Hubble 
parameter ED, 


Pi(t) = H2 




exp(-2TiJ/iLi„). 


(13) 


Here {pm)m = Pm{tin), and the time lies in the interval t € [tin,ti„ + 4/(3Ftin)]. Finally, from Friedmann’s equation 
© we hnd the energy density. 


Pit) = 




6S(e'^ — 1) + 


(Pi)i 


Hf, 


exp (-2riL/7Jin) 


(14) 


Thus we have shown how the inflation is realized, by using an inhomogeneous equation of state parameter and viscous 
coupled fluids. 

Further, we can investigate how this inflationary model conforms with recent results from the Planck satellite. We 
first calculate the slow roll parameter e, 


£ 


H 

IP 



(15) 


Acceleration takes place at e < 1. This requirement is fulfilled, when t S [tin,tin + ■\/2/(3Ftin)]- But the spectral 
parameters ns and r become negative. Consequently, the results of the recent BICEP2 experiments are not reproduced 
in the model. 
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B. Fluid model with Lo{p) = uig constant and viscosity proportional to 


In this section we will assume that the thermodynamic parameter uj(p) = wq is a constant in Eq. ([5])- We study 
the initial stage of inflation and take the function f(II) to be proportional to the square of H, 


f{H) = 

where 6 = 6/3. As in the previous section the interacting term Q has the same form as in Eq. (|5]). 
The gravitational equation of motion (jb)) takes the form: 

^-H 
3 


Wo + - ) = 0. 


The Hubble parameter is equal to 


H = 


3 (uJQ + I ) (^ - tin) - TJ- 

We find the solution of the gravitational equation of motion for matter as 

2 / \ 

35 




^o+s/e 


Pit. 


Wo - 1 + 


T^rn I - 


35 


Wo - 1 + 




(16) 


(17) 


(18) 


(19) 


where p-m = pi(tin) is the energy density of matter in the beginning of inflation. 
The dark energy density is given by the expression 


pit) - 3 ( TW + 


Wo + 1 + I 


- 


H 

W. 


^o+s/e 


Pm 


35 


Wo - 1 + 


T^in 


( 20 ) 


Thus, we have constructed an example of an inhomogeneous viscous fluid cosmology coupled with matter, applied to 
inflation. 

Further, we will use the solution (fTSl) to calculate the slow roll parameter in the inflation, 


£ = 


Wo 


6 


( 21 ) 


where wo + 5/0 > 0. As above, in order to have a regime of acceleration, one must require £ < 1, thus wq < (2/3) — 5/0. 
Another important slow roll parameter is M \ 


rj = e — 


In our case, e = rj. The power spectrum is 0 


2eH 




-£■ 


- 


^ 87r2£ ■ 

Using Eq. (1181) for the Hubble parameter, the power spectrum takes the form 

P 


K = 


3P ^Wo + 3 ^Wq + I) ~ ^in) — 7^ 

From the slow roll parameters one can calculate the spectral index Us and the tensor-to-scalar ratio r, 

rig = 1 — 6£ + 2 ? 7 , r = 16£. 

We obtain 

Ug = 1 - 6 ( Wo + ^ ) , r = 24 ( Wo + ^ V 


0 


0 


( 22 ) 

(23) 

(24) 

(25) 

(26) 


In the particular case 5 = 0, we obtain the same values for the parameters dun, (HH), and dMl) as in Ref. [43| . 

From the observations by the Planck satellite it is known that rig = 0.9603 ± 0.0073. In order to satisfy this result, 
we must require ojq + 5/0 ~ 0.00(6). Consequently, the present model can produce inflation. 
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III. QUASI-DE SITTER EXPANSION FOR INFLATION PRODUCED BY TWO COUPLED FLUIDS 


In this section we will investigate a nonviscous model for the cosmic fluid. The equation of state is taken to have 
the following inhomogeneous form [43 |. 

w(p) = -l + ai/9i/2_^^^-i/2^ ^27) 

where ai and 02 are positive dimensional constants. This model describes a quasi-de Sitter inflationary expansion. 
We will find an exact solution of this model taking into account the coupling of two different fluids. 

Let us start from the gravitational equation for energy, 

p + - 02 ) =-Q. (28) 

We take the energy exchange between energy and matter in the form 

Q = (29) 

where i5 is a positive dimensional constant, and we put the constants oi and 02 in Eq. (1271) equal to 


1 


oi = 


02 = 


Taking into account Eqs. ( 0 , dini) and dMl) we obtain the continuation equation for energy 

3 o 

= 0 . 

Ve 


The solution gives the Hubble parameter in the form 

H = 


1 


{t tin) + T7T 


(30) 

in the simple form 

(31) 

(32) 


■'.ny , 

Eurther, we find the solution for the gravitational equation of matter. 


Pi{t) = 3p'/9 e' 


2Ve 


C' + 2Ei 


' 2Ve " 


(33) 


where F,i{2'/d/H) is the integral exponential function and C an arbitrary constant, where is the integral exponential 
function and is the arbitrary constant. 

Then the expression for the energy density becomes 


p{t) = 3 




C' + 2Ei 



(34) 


Thus we have obtained a situation in which the equation-of-state parameter changes slowly during the inflationary 
period, taking into account the coupling of two fluids. 

As in the previous section, we will consider how this inflationary model confirms with the Planck observational 
data. First, let us calculate the slow roll parameters e and 77 . From the solution (l32|) we get 

e = ri = ^£, (35) 

2 VO 2 

and the power spectrum is given by the expression 


A 


2 _ 
R — 


eVo 

127r2 


The spectral index Us and the tensor-to-scalar ratio r are given by 


(36) 


rig = 1 — 3e, r = 16e. 


(37) 
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For a de Sitter expansion we may introduce the number of TV-folds, 


N = 





In our case 


N = 


2V0 




(38) 


(39) 


The inflation is viable if > 76. Therefore, the slow roll parameters change during inflation, 




and the slow roll conditions are satisfied. 
The spectral indices are given by 


ris 



r = 


16 

IV’ 


(40) 


(41) 


For N = 76 the indices are rig = 0.9605 and r = 0.221. Consequently, the results from the Planck satellite data can 
be realized from this model. According to the recent Planck data the tensor-to-scalar ration is constrained to be 
r < 0.11 (95% CL). 


IV. CONCLUSION 

In the present paper we investigated coupled-fluid cosmological models which take into account viscosity properties 
of the fluid in FRW flat space-time, in a hot universe. We considered the influence of the interaction between energy 
and matter during the inflationary very early stages of the evolution. We paid attention especially to the initial 
stage of the inflation. Studying the agreement between the theoretical inflationary models and the last results of 
the Planck satellite data, we showed that some restrictions on the thermodynamic parameters allowed satisfactory 
correspondence with the observations. 

Inhomogeneous fluid cosmology may be understood as some kind of modified gravity (for a review, see Ref. 13711) . 
owing to the fact that it may be presented as a gravitational fluid with an inhomogeneous equation of state [39j. 
Then, it follows that the unification of inflation with dark energy in modified gravity one can also achieve this 
unification in viscous cosmology, with the inclusion of dark matter. 

Our theory may be extended for the case of multiple coupled viscous fluids. The calculation may be done in the 
same way as above. 
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